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Selection problems of Z^-periodic 
entropy solutions and viscosity solutions 
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Abstract 

Z^-periodic entropy solutions of hyperbolic scalar conservation laws and Z^-periodic 
viscosity solutions of Hamilton-Jacobi equations are not unique in general. How¬ 
ever, uniqueness holds for viscous scalar conservation laws and viscous Hamilton- 
Jacobi equations. Ugo Bessi [1] investigated the convergence of approximate 
Z^-periodic solutions to an exact one in the process of the vanishing viscosity 
method, and characterized this physically natural Z^-periodic solution with the 
aid of Aubry-Mather theory. In this paper, a similar problem is considered in the 
process of the finite difference approximation under hyperbolic scaling. We present 
a selection criterion different from the one in the vanishing viscosity method, which 
exhibits difference in characteristics between the two approximation techniques. 
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1 Introduction 


We consider hyperbolic scalar conservation laws and the corresponding Hamilton-Jacobi 
equations 

( 1 . 1 ) ut +H{x,t,c + u)^ = ^, 

(1.2) vt +H{x,t,c^v^) = h{c), 


where c, h{c) G M are given constants. Here, the function H{x,t,p) is assumed to satisfy 
the following (H1)-(H4): 


(HI) H{x,t,p) : T2 X 
(H2) Hpp > 0, 

(H3) lim 

IpK-I-oo \p\ 


M, C^, 


= +CO, 
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where T := M/Z. T is identified with [0,1), and a function defined on T is regarded as 
a 1-periodic function defined on M. We say that a function f{x,t) is Z^-periodic, if it 
is 1-periodic in both x and t. It follows from (H1)-(H3) that the Legendre transform 
L{x,t,^) of H{x,t, •) is well-defined, and is given by 

L{x,t,Cl = sup{^p- 

peK 


Note that the function L satisfies 

^ M, 


(LI) L{x,t,^) : T2 X 
(L2) Lgg > 0, 

(L3) lim 

|«K+cx> lel 


= -|-oo. 


The final assumption for H is: 


(H4) There exists cc > 0 such that \Lx\ < <^{\L\ + 1). 

A flux function H satisfying (HI) and (H2) is common in continuum mechanics. A 
function H satisfying (H1)-(H4) is common in Hamiltonian dynamics, and is called a 
Tonelli Hamiltonian. Note that (H4) implies completeness of the Euler-Lagrange flow 
generated by L, and hence the Hamiltonian flow generated by H. Aubry-Mather theory 
extensively investigates the Hamiltonian dynamics and Lagrangian dynamics generated 
by a Tonelli Hamiltonian H and its Legendre transform L. Weak KAM theory unifies 
Aubry-Mather theory, Z^-periodic entropy solutions of fll.ip and Z^-periodic viscosity 
solutions of fll.2p . providing many useful tools for the analysis of the PDEs and dynamical 
systems [g], [1], [7]. 

We briefly explain the background of our selection problem in terms of the large-time 
behaviors of entropy solutions and viscosity solutions. It is well-known that the initial 
value problems 


(1.3) 


(1.4) 


Ut +H{x,t,c +u)x = ^ inTx(0,T], 

u{x,0) = u^{x) E L°°{T) on T, f u^{x)dx = 0, 

J Vt + H{x,t,c + Vx) = h{c) inTx(0,T], 

1 n(a:, 0) = v^{x) E Lip{T) on T, 


are uniquely solvable in the sense of entropy solutions and viscosity solutions, where u E 
(^^((O, T], L^(T)) and v E Lip(T x (0,T]), respectively. In the spatially one-dimensional 
case, fll.Sp and fll.4l) are equivalent in the sense that the entropy solution u or viscosity 
solution V is derived from the other if In particular, we have u = Vx (see, 

e.g., m)- From now on we assume that and that an entropy solution u E 

^^((O, T], L^(T)) means the representative element given by Vx- 

The viscosity solution v of fll.4p exists for T —)■ -|-oo and tends to a time-periodic 
viscosity solution of fll.21) with a period greater or equal to one as t ^ oo [3] (and 
the references cited there). The periodic state may depend on initial data, namely, 
time-periodic viscosity solutions with each period are not unique with respect to c in 
general. Since the entropy solution u is equal to Vx, similar large-time behaviors and 


2 















the multiplicity of periodic states hold for fll.Sp and fll.ip [2]. In particular, Z^-periodic 
entropy solutions of fll.ip and Z^-periodic viscosity solutions of fll.2p are not unique with 
respect to c in general. Note that for these large-time behaviors, h(c) must be a special 
value called the “effective Hamiltonian”. Otherwise v has linear growth in time. The 
effective Hamiltonian is in one-dimensional problems. For details on the effective 
Hamiltonian, see [2]. From now on, h{c) is assumed to be the effective Hamiltonian. 

It is common to approximate the entropy solution u of fll.Sp and viscosity solution v 
of fll.dp by the smooth solutions of the following parabolic equations with z/ > 0, the 
periodic boundary condition and the same initial data as the above, 

(1.5) < + i/(x,f,C-t-M^)a; = Z/mL, 

(1.6) + H{x, t,c + O = h^(c) + 

Within each bounded time interval, the convergence ^ u and ^ n as —)■ 0 
can be proved. This is called the vanishing viscosity method. For each initial data, the 
solutions and v'' exist for T — )■ -|-oo, and tend to time-periodic solutions of fll.51) 
and of fll.6p as f — ?■ oo, respectively. Unlike the inviscid problems, and v'^ are 
unique with respect to c and their period is exactly equal to one (to be precise, v’^ is 
unique up to constants). Note that we need to choose an appropriate constant h^(c) 
for the asymptotic behavior. h^(c) is also called the effective Hamiltonian, and tends 
to h(c) as z/ —)■ 0. From the family {u’^}u>o (resp. {v‘'}u>o (adding a constant to v^, if 
necessary)), we can take a convergent subsequence, whose limit is a Z^-periodic entropy 
solution of fll.ip (resp. Z^-periodic viscosity solution of fll.2p i [12], [1|. Since Z^-periodic 
entropy solutions and viscosity solutions are not necessarily unique with respect to c, 
an interesting problem arises: Does {h^}y>o (resp. ^.>0) accumulate on a single 1?- 
periodic entropy solution u (resp. I?-periodic viscosity solution v) as u ^ 0? If this is 
the case, how can u and v be characterized? This problem is partially solved in [12], [1]. 

Selection problems arise also in higher dimensional stationary problems 

(1.7) H{x,t,cVx) = h{c) in ¥"■, 

where the viscosity solutions are not unique in general with respect to c G M"". The 
vanishing viscosity method also works with the elliptic equation 

(1.8) H{x,t,c +v((,) = h''{c)-I in 

which is uniquely (up to constants) solvable for each z/ > 0 and c. From the family 
{l^}u>o (adding a constant to each if necessary), we can take a convergent subse¬ 
quence, whose limit is a viscosity solution of fll.7l) . The selection problem in fll.7p and 
fll.Sp is partially solved in [T] . There is another approximation method called the ergodic 
approximation with the discounted Hamilton-Jacobi equation 

(1.9) ev^-\-H{x,t,c +vD = h{c) in 

This is also uniquely solvable for each e > 0 and c. From the family {u^}e>o, we can 
take a convergent subsequence, whose limit is a viscosity solution of fll.7l) . The selection 
problem in fll.7p and fll.91) is partially solved in [11] and then almost completely in 
[6], based on weak KAM theory. Recently, another approach to this selection problem 
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is announced in [T3] based on the nonlinear adjoint method, which covers degenerate 
viscous Hamilton-Jacobi equations as well. 

The finite difference approximation is also a common technique to obtain the entropy 
solution u of fll.31) and viscosity solution v of fll.jp . which is a simple and realistic ap¬ 
proximation available on a computer. In the case of a Tonelli Hamiltonian, however, 
it is not easy to verify the stability and convergence of hnite difference schemes for 
T -|-oo and to study large-time behaviors of difference solutions for possible peri¬ 
odic states. Recently, the author announced a new approach to the Lax-Friedrichs hnite 
difference scheme based on probability theory and calculus of variations, and obtained 
new results on time-global stability, large-time behaviors, error estimates, existence of 
periodic difference solutions, etc., with lots of useful details in application of the scheme 
to weak KAM theory [18], [19]. These arguments pose a selection problem of Z^-periodic 
entropy solutions and viscosity solutions as well. A Lax-Oleinik type operator for the 
Lax-Friedrichs scheme, which is introduced in [TS], [IH], is the basic tool for the investi¬ 
gation of the selection problem. We discretize fll.jp with the Lax-Friedrichs scheme and 
fll.jp with a scheme so that the following two difference equations are equivalent: 


( 1 . 10 ) 

( 1 . 11 ) 


U 


k+l 

m+1 




(“m+^m+2) 

2 

At 

_ 2 

At 


H{x„i+2,tk,C + u’f^^2) - H{Xm,tk,C + u’ff) 
2A^ 

+ H(x^, h, c + = ftA(c), 

2Ax 


0 , 


We will give details on the discretization in Section 3. It is proved that the two schemes 
are globally stable with hxed A = {Ax, At), if A := At/Ax < Ai with an appropriate 
number Ai and if h/\{c) is chosen properly. hA{c) is also called the effective Hamiltonian 
for the difference Hamilton-Jacobi equation, and tends to h{c) as A —)■ 0. Furthermore, 
any solutions and vi+i, tend to periodic difference solutions of fll.lOp and ^m+l 
of (11.lip as the time-index k —?■ -l-oo. resoectivelv. and 'vt+i are unique (up to 
constant for ^m+l) with respect to c, and have the period one. Let ma be the step 
function derived from and let La be the linear interpolation of Then, from 
the family {'UA}Ax>o,At=AAa: (resp. {hA}Aa;>o,A 4 =AAx (adding a constant to each ua, if 
necessary)) with hyperbolic scaling-. 0 < Aq < A < Ai (A is hxed), we can take a 
convergent subsequence, whose limit is a Z^-periodic entropy solution of (11.ip (resp. 
Z^-periodic viscosity solution of (ll.2p j. Here is the selection problem: 


Selection problem. Does {'UA}Ax>o,At=AAa: (resp. {hA}Aa:>o,A 4 =AAxj accumulate on a 
single If-periodic entropy solution u (resp. -periodic viscosity solution v) as A ^ 0 
under hyperbolic scaling? If this is the case, how can u and v be characterized? 

It follows from [16] that if we take diffusive scaling, i.e.. Ax —?■ 0 with Ax^/At = 0(1), 
then ua and va tend to the unique Z^-periodic solutions of (II.Sp and (II.6p . respectively, 
and no selection problem arises. 

The purpose of this paper is to present an answer to the selection problem above, 
comparing it with the selection problem in the vanishing viscosity method [1]. The 
contributions of this work are: 


(1) This is the hrst attempt to formulate and solve the selection problem of Z^-periodic 
entropy solutions and viscosity solutions in the hnite diherence approximation. 
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(2) Mathematically and physically new aspects of the hnite difference approximation 
are made clear in contrast with the vanishing viscosity method, particularly from 
the viewpoints of scaling limit of random walks, 

(3) Our result leads to better understanding of various numerical experiments. 

2 Result 

Throughout this paper, it denotes a Z^-periodic entropy solution of fll.ip . v a Z^-periodic 
viscosity solution of fll.21) . a Z^-periodic difference solution of fll.lOl) . v^+i ^ 
periodic difference solution of fll.lip . ma the step interpolation of and the linear 
interpolation of ^m+l- In order to specify the value c, the notation , v^\ 

Um(c), h^+i(c), etc., is sometimes used. Dehne quotient maps 

pr : M 9 a; X mod 1 G T, 

pr : 9 (x, s) i-G- (x mod 1, s mod 1) G T^, 

pr : 3 (x, s, y) ^3 (x mod 1, s mod 1, p) G x M. 

In this paper, we use the term “projection” for the operation of pr (it does not mean 

“T X R 9 (x,p) !-)■ X G T”). 

According to weak KAM theory, we have a characteristic curve 7 : (—cxo,t] —)■ R of 

p(c) Qj. p(c) _ yh) 

7 (t) = X and 'y'{s) = Hp{'y{s), s, c + s)), s < t 

for each x, f G R. Furthermore, 7 falls into a global characteristic curve 7 * : R —)■ R as 
s -3 —00 (otherwise 7 is a part of a global characteristic curve). The curve 7 * has the 
rotation number 

lim = h'{c). 

|s|^oo S 

For each c, the set 

;= lj{pr(7*(s),s) I s G R} C 
7* 

is called the Aubry-Mather set for c, where the union is taken over all the global char¬ 
acteristic curves 7 * of all Z^-periodic entropy solutions or viscosity solutions with c. It 
is known that, if the rotation number h'{c) is irrational, and are unique (up to 
constants for [7]. Hence, we consider the case where h'{c) is rational. In particular, 
we deal with the situation where the global characteristic curves yielding form 

smooth hyperbolic stable/unstable manifolds, and the graph of each 

graph(c -I- := {(x, t, c -1- -u(‘^)(x, t)) | x, t G T}, 

consists of parts of these stable/unstable manifolds. It is easy to see an example of such 
a situation with non-uniqueness of it and v through explicitly solvable problems given 
by Hamiltonians of the form H{x,t,p) = |p^ — F{x). Here are the assumptions for our 
selection problem: Throughout this paper, A = At/Ax is hxed with 0 < Aq < A < Ai, 
where Ai is from the stability analysis of (ll.lOh and (ll.llh and Aq is arbitrary. 
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(Al) For all c G [co,Ci], the Aubry-Mather set consists of common hyperbolic 

periodic orbits jI, ... with a rational rotation number p/q, q E N,p E Z, where 
0 < 7 i( 0 ) < 72 ( 0 ) < • • • < 7 /( 0 ) < 1 and {pr{'-y*{s), s) | s G M} 7 ^ {pi^hjis), s) \ s E 
R} fori ^ j. 

(A2) If c = Co (resp. C\), graph{c+vf'A'^ coincides with the lower separatrix (resp. the up¬ 
per separatrix) formed by the projected -stable/unstable manifolds of ■jl,... 

(A3) LetU)^^{5) denote the unstable manifold 0 / 7 * restricted to 

U\5) ■.= {{x,t)eM^\\x--f:{t)\<6}. 

Let C{x,t) be a C"^-generating function of the local unstable manifold , i.e., 

{{x,t,Cx{x,t)) I {x,t) e U^{S)} = Ul^{5), and set 

r*:= r{CUA:{s),s)-X\U^:{s),s)}ds, t 
Jo 

Then, mirij F* is attained by only one i = i*. 

(Af) If the limit of a convergent subsequence satisfies graphic + D 

prlAl^^{5) for small 5 > 0 , then is equi-Lipschitz with respect to x on U'‘{5), 
i.e., 

l^m +2 (^) ut{c)\ < 60 ■ 2 Aa; for all {xm,tk), {xm+2,tk) e 
with Oq independent of A, c and the choice of the subsequence. 


We refer to known results on and in Section 3. Before stating our result, we 
recall the Peierls barriers, which play an important role in Aubry-Mather theory. For 
each c, the Peierls barrier h^\x, t] y, r) : x ^ R is dehned as 


h^pKx,t-,y,r) 


lim inf 

TGN,r->oo 


fr+T 


inf 


{L(k(s),s, k'(s)) 


ck'{s) -\- h{c)}ds , 


where the inhmum is taken over all absolutely continuous curves k : [t, r -|- T] —)■ T with 
k{t T) = y and K{t) = x. For details on the Peierls barrier, see [13], [9]. Our main 
result is: 


Theorem 2.1. Suppose that (Al)-(A4) hold and that c E (co,Ci). Add a constant to 
so that h^^( 7 **( 0 ), 0) = 0. Then, as A —)■ 0 with 0 < Aq < At/Ax = A < Ai, 


(1) converges to hp'^^( 7 **( 0 ), 0 ; •, •) uniformly, 

(2) converges to ■, ■))x pointwise a.e. and in the C^{T] L^(T))- 

norm. A geometrical characterization of u)^'> is that there exists 5 > 0 such that 
graphic -|- contains prU/^^{5) only for i = i*. 


We compare our result with the one in the vanishing viscosity method [1|. Our setting 
given by (Al) and (A2) is essentially the same as the one in [1|. In both of the cases, 
the basic ingredients are stochastic Lax-Oleinik type operators for fll.bp and fll.lll) . and 
estimates for the rate of the law of large numbers realized in the operators as z/ ^ 0 and 
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A —)■ 0. In |1], the standard framework of the stochastic Lax-Oleinik type operator for 
viscous Hamilton-Jacobi equations HD] is used. We will see our stochastic Lax-Oleinik 
type operator for fll.lip in Section 3. Although the mechanisms of the selection are 
similar to each other, the selection criteria are different. Roughly speaking, our case also 
relies on the so-called numerical viscosity of the Lax-Friedrichs scheme, which causes 
diffusive effect similar to that of the artihcial viscosity term However, the speed of 
the diffusion caused by the numerical viscosity is Ax/At = and due to hyperbolic 
scaling, it is hnite for A —)■ 0. This leads to the different selection criterion. In fact, 
(A3), which gives our selection criterion, is different from the one in |1]. For the vanishing 
viscosity method, the values 

f’:= 

Jo 

appear, and h^(c) selects hjf^( 7 -*,( 0 ), 0; •, •) with i* that minimizes f*. In our case, the 
selection criterion contains the discretization parameter A(> Aq > 0) and we may change 
the selection by varying A, which means that our hnite difference method with hyperbolic 
scaling is more “hexible” or “closer” to the exact hyperbolic problem than the vanishing 
viscosity method. Our method does not always select the physically natural solutions 
selected by the vanishing viscosity method. If we choose A > 0 small enough, then our 
selection is the same as the one in the vanishing viscosity method. This is an interesting 
difference between the two methods, since they are often regarded as mathematically 
similar due to their diffusive effects. At the current stage, the regularity assumption 
(A4) plays an essential role yielding an appropriate estimate for the rate of the law 
of large numbers, whereas there is no such assumption in [1]. This is due to different 
structures of the stochastic Lax-Oleinik type operators. In the case of the vanishing 
viscosity method, it is enough to deal with the standard Brownian motions. In our case, 
we have to investigate continuous limits of space-time inhomogeneous random walks with 
hyperbolic scaling. It is not easy to prove the law of large numbers and to estimate its 
rate for such random walks im. As far as the author knows, justihcation of (A4) is an 
open question, though computer simulations imply that it would be true. Particularly, 
in the case of entropy solutions with hnite number of shocks, their Lipschitzian parts 
seem to be approximated in equi-Lipschitzian ways (see, e.g., m)- 


3 Preliminaries 

We state several known facts on Z^-periodic entropy solutions, Z^-periodic viscosity 
solutions and our difference schemes, as well as space-time inhomogeneous random walks 
arising in the difference equations. 

3.1 Z^-periodic viscosity solution 

Let V = be a Z^-periodic viscosity solution of fll. 2 p which is periodically extended to 
Then, v satishes the following equality for each x, f G M and t < t (the deterministic 
Lax-Oleinik type operator): 

(3.1) h(a:,t) = inf | /" ^^'^^( 7 ( 5 ), s, 7 '(s))ds + h( 7 (r), r) 1 + h(c)(f - r), 

'yeAC,'rit)=x J 
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where AC is the family of all absolutely continuous curves 7 : [r, t] ^ M and 


■= L{xA: 0 - ci 

is the Legendre transform of H (x, t, c + •). There is a minimizing curve 7 *, which is a C^- 
backward characteristic curve of h, and solves the Euler-Lagrange equation generated by 
LA'> . Moreover, v is differentiable with respect to x on the minimizing curve, satisfying 
the relation 

(3.2) 7 *'(s) = iLp( 7 *(s),s,c + ha;( 7 *(s),s)) for r < s < t. 

Note that, if Vx{x,t) exists, s = t is included in fl3.2l) and 7 * is the unique minimizing 
curve for fl3.ip . Each Z^-periodic viscosity solution belongs to Lip(T^), and each 
Z^-periodic entropy solution u^) belongs to Lip(T; L^(T)). We have and 

(3.3) v^^\x,t) = U^‘^\x,t) — [ [ {V'j;^\x,T) + H{x,t,c +U^‘^\x,T))}dxdT 

Jo Jo 

+h{c)t, 

where U^^\x,t) ■= Jq A)du- For more details, see 0 , 0 , 0 , 0 - 


3.2 Discretization 


Here are some details of our discretization and results on the difference equations ob¬ 
tained in [18], [19]. Let K,N be natural numbers with N < K. The mesh size 
A = (Ax, At) is defined by Ax := {2N)~^ and At := {2K)~^. Set A := At/Ax, 
Xm '■= mAx for m G Z and t^ := kAt for A; = 0,1, 2,.... Let (AxZ) x (AtZ>o) be the 
set of all {xm,tk), and let 


Qeven C (AxZ) X (AtZ>o) (resp. Qodd c (AxZ) x (AtZ>o)) 


be the set of all (x^, t^) with A; = 0,1, 2,... and m G Z such that m -|- A; is even (resp. 
odd), which is called the even grid (resp. odd grid). For x G M and t > 0, the notation 
m(x), k{t) denotes the integers m, k for which x G [xmiXm + 2Ax) on Qeven or Qodd and 
t G [tfc,tfc -|- At), respectively. Note that m(x) on Qeven and m(x) on Qodd are different 
for the same x. We discretize fll.ll) on Qeven by the Lax-Friedrichs scheme as fll.lOp . We 
also discretize (USD in Qodd as fll.lip . We introduce the following notation: 


Dtw 


k+l _ 

m • 


W. 


k-\-l 


w 


k 

m — 1 


-\-W 


k 

m+1 


2 


At 


DxW 


k 

m+1 * 


wt+i 


— W 


k 

ra—1 


2 Ax 


Note that fll.lOp and fll.lip are equivalent. In particular, if + satisfies fll.lip . then 
< := Dxvt+i satishes fll.lOp . Throughout this paper, we follow the convention that 
the sum of the super-script and sub-script of variables defined on Qeven (resp. Qodd) is 
always even (resp. odd) in the notation h^+i, etc. We say that and 

k>m+i nre Z^-periodic, if they satisfy 


^.k — ^.k-\-2K — ^.k 

'^m±2N ~ “m ~ “m 


^,k _ ;r.k+2K _ -k 

^m+l±2N ~ ^m+l ~ ^m+1 


for all m, k. 






















Let ma be the step function derived from a Z^-periodic difference solution i.e., 

UA{x,t) := for {x,t) G [Xm-l,Xm+l) X [4,4+i)- 

Let va be the linear interpolation with respect to the space variable derived from a 
Z^-periodic difference solution i-e., 

va{x, t) := + -Dxh^+i • {x - Xm-i) for (x, t) e a^m+i) x [4,4+i). 

Note that va{x, •) is a step function for each hxed x and that (ua)* = ua- 

We sometimes use a phrase like “a convergent subsequence {ua} which tends to v as 
A —)■ 0”, meaning “a convergent sequence {hA^jjGN with —)■ 0 which tends to v as 
3 oo”. 

Proposition 3.1 f[TU]). Let [co, Cl] he an arbitrary interval. There exist Ai > 0 and 
(5i > 0 such that for any A = (Ax, At) with Ax < (5i and X = At/Ax < Ai we have the 
following statements for each c G [cq, Cij: 


(1) There exists one and only one number hA^c) for which M.ll\) has the unique (up 
to constants) I?-periodic difference solution vt+i- 

(2) There exists the unique Tf-periodic difference solution = D^v/n+i of U.10\) . 
which is uniformly bounded with respect to m, k, A and c with the entropy condition 


ut+2 - uj 

2Ax 


< M 


for all k, m, 


where M > 0 is independent of A and c. In particular, the stability condition 
(CFL-condition) is verified: 


\Hp{Xm,tk,C + U^)\ < X^^ < X h 
If c < c, then c + < c + u^{c) for all k, m. 


A —)■ 0 under hyperbolic scaling, i.e., Q < Xq < X < Xi, 

(3) hA^c) converges to the exact effective Hamiltonian h{c) uniformly on [co,Ci] with 
the order y/Ax. 

(4) Adding a constant to each va if necessary, we can subtract a convergent subse¬ 
quence, which tends to a function v uniformly, v is a 1?-periodic viscosity solution 

ofn- 

(5) Let {ua} be the convergent subsequence in (4). Then ua = {va)x converges to 
the 1?-periodic entropy solution u = Vx pointwise except for the points of non¬ 
differentiability ofv and also in the ^^(T; L^{T))-norm. In particular, ua converges 
to u uniformly in the outside of an arbitrary neighborhood of shocks. 


We remark that u and v are obtained as a result of large-time behaviors of any other 
difference solutions [T9] . 
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We show the stochastic and variational structure of fll.lll) . First, we introduce space- 
time inhomogeneous random walks in Qodd^ which play the role of “characteristic curves” 
for fll.lOl) and fll.lll) . For each point G Godd, we consider backward random 

walks 7 that start from Xn at ti+i and move by ±Aa; in each backward time step At: 

7 = .,«+!, = Xn, 7 ^^^ - 7 ^ = 

More precisely, for each G Qodd we introduce the following objects: 

X . I ^ Goddi |^m+l ^n\ — H“ 1 /u)Ax} for ^ 1 , 

G := IJ (A" X {4}) C Qodd, 

I' <k<l-\-l 

^ : G 3 {xm+i,tk) ^ ^ A = At/Ax, 

p : G 3 (xm+1,4) Pm+i ■= 2~ ^ 

p : G 3 (xm+1,4) '-t Pm+i 2 ^ [O’ 1]’ 

7 : {/', I' + 1,I + 1} 3 k i-3 'y^ e X^, 7 ^+^ = x„, 7 ^+^ — 7 ^ = ±Ax, 
f 2 : the family of the above 7 . 

The value p^+i (resp. p^+i) is regarded as the probability of transition from {xm+i,tk) 
to (xm+i + Ax,tk — At) (resp. from {xm+i,tk) to (x^+i — Ax,tk — At)). The function 
^ is a control for random walks, which plays the role of a velocity held on the grid. We 
dehne the density of each path 7 G fl as 

h(7) := n F(A^7^■^), 

I' 

where p( 7 ^, 7 ^“^) = ^^( 7 '=) Pm( 7 '=)) — 7 ^“^ = —Ax (resp. Ax). The density 

/i(-) = p (-;0 yields a probability measure for 11 , i.e., 

prob{A) = /i(7; 0 A C H. 

7 GA 

The expectation with respect to this probability measure is denoted by namely, 

for a random variable / : H —)■ M we have 

^M';o[/( 7)] = 5^h(7;0/(7)- 

7 G 0 

We use 7 as the symbol for random walks or a sample path. If necessary, we write 
7 = 7 (xn,fi+i;C in order to specify its initial point and control. 

We state several important results on the scaling limit of such inhomogeneous random 
walks, obtained in im. Note that the variance := i?^(.;^)[| 7 ^ — 7 ^P] with 7 ^ := 
F'^(.;^)[ 7 ^] is of the order 0 ( 1 ) for inhomogeneous random walks in general, whereas 
it is of the order 0(Ax) for space homogeneous cases (see Remak 3.3 of [E]). Let 
7 ( 7 ) = {v^h)}k=i', //+!,,7 G H be a random variable that is induced by a random 
walk 7 = 7(xn,tz+i;0 as 

(3.4) r^^+C= 7 '+\ y( 7 );= 7 ^+i- ^ ,tk>)At for/'< </. 

Set := i?^(.;^)[| 7 ^ — 7 ^( 7 )P] and := i?^(.;^)[| 7 ^ — 7 ^( 7 )]] ioi I' < k < I + 1. 
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Proposition 3.2 ([IZ]). (1) For any control^, we have 


<a^ < 

A 


(2) If a control f satisfies the Lipschitz condition at 7 ; 


with C ■= Cm + 


t/c 

m(7^) f—,k 


m+1 


^*1 < 0\Xm+l - 7 


;r,fc| 


2 Ax (1 ~ ^m( 7 '=)) for dll k,m, then we have 

46(ii+i—4) 


k ^ 0 

a'" < — 


40A 


-Ax. 


If we take the hyperbolic scaling limit A = (Ax, At) —)■ 0 with 0 < Aq < A = At/Ax < 
Ai, then and always tend to zero with 0(Ax), and so does with the above 

Lipschitz condition. (A4) will be used to verify the Lipschitz condition. 

Let C+i be a Z^-periodic difference solution of fll.lip that is periodically extended 
on the whole of Qodd- Then C+i satishes the following equality for each n, I and I' < I 
(our stochastic Lax-Oleinik type operator) [18]: 


(3.5) 7+' 


illf £■„(.;£) 


5^ 

I' <k<lY-l 


+ h/\{c)(ti^i — tifi, 


where the inhmum is taken over all controls bounded by A We can hnd the unique 
minimizing control This satishes with Ai in Proposition 13.11 

let+'|<Ar'<A-' and CC = c + D^C+i)- 

The equality (13.Ah is the key tool for our proof of Theorem 12.11 


3.3 Construction of Z^-periodic solution 

We observe how to construct Z^-periodic entropy solutions and viscosity solutions with 
the given position and number of singularities under the assumptions (Al) and (A2). 
This is necessary in Section 4. For k = 0,... ,q — 1, dehne 

:= prbiik),C+ii^)], i = l,...,/-l, 

:= P?^[7/(^),7i(^ +1)], if ? > 1, 

0 /,o := pr[ 7 ;( 0 ), 7 i( 0 ) + 1], if g = 1. 

Note that = T. Let (resp. C~i/) be the segment of the upper (resp. 

lower) separatrix restricted to 0*^ x {t = 0}. For each i, choose arbitrarily either 
or {Cj Q,..., Cj g_^}. Then, our choice yields a curve C = {(x, 17o(x)) | x G 
T} that consists of or on each Qi^k- Let c := Uo(x)dx and let Uq := Uq — c. 
Then, the solution of fll.3p with this c and Uq is a Z^-periodic entropy solution without 
any shock. 

Suppose that ..., Cfiq_fi\ was chosen in the above C for some i. Fix arbitrarily 

/co e {0,..., g — 1} and take go ^ Let C = {(x, Uq{x)) | x G T} be the curve such 

that C = C on T \ Qi^ko and C switches from Cfi^ko Vo as x increases in 0j,fcp. 
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Let c := /g ilo {x)dx and let uq := Uq — c, where Ac := c — c —)■ 0 as yo tends to the right 
boundary of Qi^ko- Then the solution u of fll.3j) |.,g-.-g is a g-periodic entropy solution 
with a single shock. In fact, u{-, 1) has only one jump at yi G u{-,2) has only 

one jump at 1/2 G ©*,^ 0 + 2 ; • • • ; u{-, q) has only one jump at yq G ©i,fco+<? = Vq = Vo 

due to the conservation law {c + u{x, t)}dx = c for all t. 

With the above notation, consider the curve C = {{x,Uq{x))\x G T} such that 
(7 = (7 on T \ Uo<fc<q©j,fc and C switches from to at yk as x increases in 

©i,fco+fc for A; = 0,..., g — 1, where ko + k is replaced by fco + A; — g, if feg + A; > g. Let 
c := Jq Uo{x)dx = c — qAc and let ho := Ifo ~ c- Then, the solution u of fll.3p |.,p-.-p 
is a 1-periodic entropy solution, i.e., a Z^-periodic entropy solution. 

In this way, we can construct a Z^-periodic entropy solution for each c, choosing the 
position and number of shocks, and therefore the corresponding viscosity solution via 

(O. 


4 Proof of result 


First, we briefly state our strategy. The key point is that our discretization scheme is 
of the hrst order, namely, for a C^-solution v of fll.2p such that \Hp{z,t,c + Vx{x,t))\ < 
Aj"^ < A“^ and for := v{xm+utk) dehned on Qodd, we have 

/\ 

T H{Xmi ^ki C DxV^j^-^ ^(c) Afc) “ 1 “ six^n^ tk)i 

A{x, t) = Vxxix, t) - X^Vtt{x, t), e{Xjn, tk) = o(Ax). 

In fact, this follows from the following Taylor expansions around {xm,tk) 

/\ rp 

= Vt{xm, tk) - —A{xm, tk) + o(Aa:), = Vx{xm, tk) + o(Aa;). 

It follows from a similar reasoning to obtain (j3.5p (see Proposition 2.2 of [18]) that 
satishes for any n, I, /' < I, 


(4.1) = infBM-4) 


Ax 




l'<k<l+l 


+£(7^,4-i)}Af+ + h{c){ti+i - ti>). 


By comparing fl3.5p and fl4.ip near 7 *, we obtain upper and lower estimates for h{c) — 
h/\{c) in terms of A and so on. Since the law of large numbers holds for the minimizing 
random walk, where it tends to 7 *, we obtain P* through A in the estimates. This leads 
to our criterion. Of course Z^-periodic viscosity solutions v are only Lipschitz, and our 
argument is more complicated with additional terms. In the estimates for h{c) — h/\{c), 
the term AxP* appears as well as other terms. In order to keep AxP* as the principle term 
in the estimates, we need sharper estimates for the others with the aid of Proposition (3^ 
The main difhculty is that Chebychev’s inequality in our random walks yields estimates 
of the oder 0(Ax) at best ~ the same order as AxFk We will refer to this point in the 
proof. 
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We say that a Z^-periodic entropy solution of or Z^-periodic viscosity solution 
of fll.2p has “transition points” at 7*, if graph(c + or graph(c + contains 
prUloci^) for small <5 > 0 (“transition” means that the graph continuously switches from 
the lower separatrix to the upper one along 7*). The following fact is used in |1] , but for 
the reader’s convenience we repeat it with a proof. 

Proposition 4.1. Let c G (co,ci), i G {1,2 ,and he a l?-periodic viscosity 
solution of U.^) with transition points only at 7*. Then, such a is unique (up to 
constants) and coincides with h^\'y*{0), 0; •, •). 

Proof. Let v = v^‘^\v = be Z^-periodic viscosity solutions of fll.2p with transition 
points only at 7*. Adding a constant if necessary, we have 

(4.2) ^(^;(0),0) = h(7;(0),0) = 0. 

We will prove v = v. Take any {y, r) G T^. Since v has transition points only at 7* and 
every minimizing curve yields a trajectory on the projected stable/unstable manifolds, 
we have the following two cases and in both of the cases we obtain v{y, r) — v{y, r) > 0: 
Let a G Z and let 70 : (—00, r] —)■ M be a minimizing curve for v{y, r). 

(i) {pr’(7o(s), s, c + h3;(7o(s), s))}s<t- is on the projected unstable manifold of 7*. Then, 
there exists 6 G Z, \b\ < q for which 70 satishes 

pr-yQ{aq + 6) —)■ 7*(0) as a ^ — 00 . 

Hence, it follows from fl3.ip that v{y, T)—v{y, r) > h(7o(ag+6), aq+b) —v^'jQ^aq+b), aq+b) 
for any a. It follows from fl4.2p that h(7o(ag + b), aq + b) — '0(7o(ag + fe), ag + 6) —)■ 0 as 
a — 00 . Therefore, we obtain v{y, r) — v{y, r) > 0. 

(ii) Otherwise. Then there exist Xi,... ,Xk G {pr'y*{t) \ t = 0,... ,q — 1, j ^ i}, k>l 
for which we have the heteroclinic chain connecting y and 7*(0) that consists of the 
following curves: 

70 : (— cxo, r] M, the minimizer for v{y, r), pr'yo{aq) —)■ Xi as a ^ — 00 ; 

71 : M —)■ M, prjii^aq) —^ X 2 (resp. Xi) as a ^ —00 (resp. +00); 

(4.3) ! 

7fc_i : R —)■ M, pr7fc_i(ag) —)■ Xk (resp. Xk-i) as a —)■ —00 (resp. +cxo); 

7fc : R —)■ R, pr'jkiaq + 6) —)■ 7*(0) as a ^ —00 with 6 G Z, \b\ < g, 

pr7fc(ag) ^ as a —?■ +00, 

with 

7'(s) = Hp{-fj{s),t, c + Vxi^jis), s)) for j = 0 ,..., k, 
which means that restrictions of 70 ,..., 7^ within any hnite time interval are minimizing 
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curves for v. Hence, we have 


v{y,T) = / {L(^)( 7 o(s),s, 7 o(s))+ /i(c)}cis + u( 7 o(ag),ag) 

J aq 

= lim [ {L('=)(7o(s),s,7o(s))+ /i(c)}hs + h(xi,0), 




u(a:i,0) = 


lim u(7i(ag),ag) 

a—>-+oo 


= lim 

a—>-+oo 

= lim 

a—>-+oo 


faq 


.J —aq 
raq 


{L(^)(7i(s), s, 7^(5)) + h{c)}ds + u(7i(-ag), -aq) 
{L(‘=)(7i(s), s, 7j(s)) + h{c)}ds + v{x2, 0 ), 


—aq 


v{xk, 0 ) = lim v{'yk{aq),aq) 

a—>-+oo 


= lim 

a—>- + 00 


raq 


{L('')(7fc(s), s, 7fc(s)) + h{c)}ds + v{-fk{-aq + b), -aq + b) 


_J —aq-\-b 

/ aq 

{L('^)(7fc(s),s,7^(s)) + h(c)}cis + h(7*(0),0). 

-aq+b 

Therefore, we obtain 

( 4 . 4 ) v{y,T) = lim [ {L^^\jo{s), s^iois)) + h{c)}ds 


' aq 


k-l 


/‘“'J 

+ ^ lim / {L('^)(7 j(s),s,7'(s))+ h(c)}hs 


a—>-+oo , 
j=i ^-<^1 

/ aq 

{L('=)(7fc(s),s,7fc(s)) + h(c)}hs + h(7;(0),0). 

-aq+b 

On the other hand, the variational property of v implies 

v{y,T) < lim [ {L('^)(7o(s),s,7o(s)) + h(c)}ds 


' aq 


k-l 


i-aq 

+ Y^ lim / {L('=)(7 j(s),s,7'(s)) + h(c)}ds 

/ aq 

{L('')(7fc(s), s, 7fc(s)) + h{c)}ds + h(7**(0), 0). 

-aq+b 

Thus, we conclude that v{y,T) — v{y,T) > h(7*(0),0) — '&(7*(0),0) = 0 . 

Since v has transition points only at 7*, we have the same argument, obtaining v{y, r) 
viy, t) < 0 and v = v. 

Next, we see that v{y,T) = h^\'y*{ 0 ), 0 -,y,T). Note that 


( 4 . 5 ) hj,‘=)(7;(0), 0 ; y, r) = lim inf 

^ TgN,T-)-oo 


inf y {l1‘=1(7(s),s,7'(s)) + h(c)}cis 
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where the inhmum is taken over all absolutely continuous curves 7 : [—T, r] ^ M with 
7('^) = y cind pr'y{—T) = 7*(0). Let 7^ : [—T,r] — ?■ M be a minimizing curve for the 
inhmum in fl4.5p . Then, we have for any T G N, 

v{y,T) < j {L(^)(7r(s),s,7r(s))+ h(c)}ds + h(7;(0),0) 

= j {L^''\' 1 t{s),s,'Jt{s)) + h{c)}ds. 

Therefore, we obtain v{y,T) < hp'^^(7*(0), 0; ?/, r). 

Dehne 7 : [—{2k + 1)0^ + 6, r] —?■ M with a G N, £ > 0, 70 ,..., 7^ and b in fl4.3p as 


7 o(s) for s G Jo = [-aq + e, r], 

linear line connecting 70 (—aq' + s) and 71 (aq — e) 

for s G Jo = [—aq — e, —aq + e], 

7 i(s + 2aq) for s G Ji = [—3aq + e, —aq — e], 

■* linear line connecting 71 (—aq + e) and 72 (aq — e), 

for s G Ji = [—3aq — e, —3aq + e], 

7 fc(s + 2kaq) for s E Jk = [—{2k + l)aq + b + e, —{2k — l)aq — e], 

linear line connecting 7 fe(—aq + b + e) and 7 *(—aq) 

for s E Jk = [—{2k + l)aq + b, —{2k + l)aq + b + e]. 


Then, we have 


^?H7r(0),0;|/,r) < hminf 


aGN,a->-+oo J_(^2k+l)aq+b 


{L^''\'y{s), s, 7 '(s)) + h{c)}ds 


aEN,a^+oo 
k 


/ {^^"H 7 (s),s, 7 '(s))+ h(c)}Js 
j=o 7 


+ Y1 / {^^''n7(s),s,7'(s)) +/i(c)}Js 
j=o JJj 


We can take e —)■ 0+ according to a —)■ +cxo so that the term on the third line tends to 
zero. The term on the second line tends to the right-hand side of fl4.4l) as e —)■ 0-I-. Thus, 
we conclude that hp'^^(7*(0), 0; ?/, r) < v{y^ r), obtaining v{y^ t) = h^\'y*{0), 0 ; y, r). □ 


The following proposition is a key fact: 

Proposition 4.2. Let c E (co,Ci). If the limit of a convergent subsequence 
has transition points at 7 * and 7 *, then T* = P . 

Proof. Set V = and wa = ■ There exists <5 > 0 for which we have graph(c-l- Vx) P> 

prUl^c{^),prUloc{d) and v is smooth on U'^{5)UU^{5) including its boundary. Throughout 
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this paper, ei(A), e 2 (A),... are numbers that tend to 0 as A —)■ 0 , and 61 , 62 ,... are 
constants. Since {va)x Vx a-e. as A —)■ 0 , we have with small e > 0 , 

(4.6) sup |hA — h| < ei(A) and sup |(hA)x — < ei(A). 

IR2 U'{ 5 +e)UU^ 5 +e) 

Now we investigate the solutions around 7 *. The same investigation is possible around 
7 *. Take a Z^-periodic C^-function n : —)■ M with the derivatives uniformly bounded 
in such that 


v\uiiS) = v\ui{s) and \Hp{x,t,c +v^{x,t))\ < 


Set 


g{x,t) := Vt + H{x,t,c + v^{x,t)) - h{c), 
where g = 0 on W{6). Then, •= v{xm+i,tk) solves the difference equation 

(4.7) DtV^^ + H{Xm,tk,C + = h{c) - —A{Xm,tk) + giXmAk) +e{Xm,tk), 

A{x,t) = Vxx{x,t) - X\t{,x,t), e{xm,tk) = o{Ax), 


and therefore it satishes the equality for any < I 


(4.8) u 


Z+1 


= infE^(,^) 


Ax 




I' <k<l+l 


+ 5'(7 Ak-i) + £{,1 Ak-i)}At + v^^^ii-^ + 6 .(c)(b+i — t//). 
Let denote the difference solution that yields ua. Define tCm+i oa Qodd as 


r h^+i onU\5)-.= U\5)AGodd 

:= \ ^m+l + ^m*(fc)+l " <-(fc)+l ^r Ul* {k) < 171, 

( vi+i + - < 7 fc)+i for m < m,{k), 


where {xm*{k)+i,tk) and (xm,(fc)+i,t^) stand for the right end and left end of W{6)\t=t^, 
respectively. Note that difference between space variables of the nearest two end points 
of f/*( 6 ) is ±Ax, due to | 7 *'(s)| = |iLp( 7 *(s), s, c + h 3 ;( 7 *(s), s))| < A“^. It follows from 
(14. 6 p that 

(4.9) |w^+i - n^+il < ei(A) on g^dd- 


Set 


fi^XmAk') . HAmytky C 6,a(c), 

where f{xmy tk) = 0 for m such that {xm±iytk) G D*( 6 ). Then, tc^+i solves the difference 
equation 


(4.10) Dtw’^^ + H{xmy tky c + J = hA(c) + f{x^, tk) 

and therefore it satisfies for any < I 


(4.11) w: 


Z+1 


= infE 


A-X) 


Y, + /( 7 ‘. 4 -.)}Ai + 


I' <k<l-\-l 

+hA(c)(tz+i - V). 
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We will show 


(4.12) 


- g(x„,tt)\ < e2{&) ong even' 


For denoted by m*, we have the following two cases: 

(i) (xm*,4+i) ^ U''{6), which implies (a;m.+2,4+i) e Then, for m < m*, we 

have 


^xU!m+l 


DtW, 


2Ax 


{{v^+i + - {vt_i + - v^.+i)} 


— D 


k-\-l 

m 


1 

At 


— _ <; + l _L + l 


v: 


k+1 

' m ,+2 


I ^r) 


+ Xm. + l ~ ^t. + l + ^t+1 X- ..m. + l ^ + 

= B,«y' + Ae*tt2 - A«St2 + A-'lIl.C+a - 


Hence, with fl4.7p and fl4.10p . we obtain 


(4.13) f{Xjn^tk) diXmitk) -^(^m*+2) ^fc) C T Fr(Xm*+2) ^fc) C + 

/\ '7^ 

+ A (Zl3;h^^_|_3 — tk) + £(Xm, tfc) 


2 A 


Ax 


T 2^ ^(^m»+2) ^fc) ^{Xmt+2itk)- 


(ii) (xm.,tfc+i) e Then, for m < m*, we have 


DxWm+l 






Au-y 




— {vt-i + v^,+i - - t^™,+i)} 


= D,<+‘+Aiy - Aiy+ 



Hence, with (14.7p and (I4.10p . we obtain 

(4.14) f {Xm^ ^k) di^my tk) H{Xmt ,) C + Dx'Vxa^^]) H{Xmt ,) C + Dx'Vxa,+]) 

{^xVmt,+l ~ T ~ ^{^my tk) 

2^ ) ^fc) "F ) ^fc) T ytk)y 

where g{xm,ytk) g{xm,+ 2 ytk) = 0 as A -)■ 0. 

Similar calculation is possible with m*{k). Therefore, noting that DxV^j^^ = Vx{xmy 4) + 
o(Ax) and T>xh^+i = {.Vj\)x{xmytk)y we obtain fl4.12p through fl4.6p . 

We state the law of large numbers for the minimizing random walks for (14.81) and 
(14.lip , where they tend to genuine minimizing curves for exact viscosity solutions. 


17 



















Lemma 4.3. (1) For each x & X := {x e M | |a: — 7*(0)| < 5 < 5}, t < 0 and e > 0, 
define 

C-{r) := {L^^\r{s), s, r'(s)) + g{r{s), s)}ds + v{r{t),t) + h{c){-t), 

7*; minimizer for inf >^(7), 

'y£AC,"f(0)=x 

Dx ■= : [t, 0] —)■ M I |r'(s)| < A~\ r(0) = x, C{r) — £(7*) < e}. 

Then, we have 

sup sup II r — 7* |lco([f,o])~^ 0 as e —)■ 0. 

xGX 

(2) Let rj^i^j) denote the linear interpolation of For each Xn+i G X such that 

(x„+i,0) G Qodd, v' > 0 and the minimizing random walk 7 = 7(x„+i,0;^*) within [tz',0] 
define 

:= {7 e I II r/A(7) - 7x„+i l|co([v,o])< 

Then, we have for each fixed v > t], 

inf proh{Vt\{xn+i)) as A —)■ 0. 

^ri+1 

In particular, the linearly interpolated averaged path 7^ of the minimizing random walk 
f7A is also the average ofg^{'j). See Theorem 3.2, m) satisfies || 7 a- 7 ^„+i ||co([v,o])^ 
0 as A —)■ 0, uniformly with respect to Xn+i- 

Remark. Since f(x, 0) = inf^g^c',7(o)=x-^(7) and v = v, g = 0 on U\6) with dif¬ 
ferentiability of V, the minimizer'-^* for m.i.y^Ac,'y(o)=x coincides with the unique 
minimizer^* for v{x,0) within a small time interval [e, 0]. Since the dynamics on the 
unstable manifold implies that 7*(s) tends to 7*(s) as s ^ — 00 , the points (7*(s),s) 
never come out ofUfid). Hence 7* is equal to 7* on [t, 0], and is the unique minimizer 
for AC,'y(o)=x £^{ 1 ) ■ It follows from hyperbolicity of 'j* that there exists bi > 0 for 

which we have for s < 0, 

(4.15) |7:(s) - 7.*(ii)l < hyo) - 7.'(0)|e''‘'. 


Proof. (1) Proceeding by a proof via contradiction, we assume that there exists a se¬ 
quence Sj —)■ 0 as j —)■ cxo for which we have 62 > 0 such that 

sup sup II r — 7* llcoA 362 for all j. 

We can take Xj G X and G such that 

Iki - llco> 262 for all j. 

We have a convergent subsequence of {xj}, still denoted by {xj}, that tends to Xjj G A 
as j -A- 00. Since 7*/(0) = Hp{xj,Q,c F Vx{xj,Qi)) 7^/(0) = i4p(xj, 0, c-1 -0)) 
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as j oo, we have 7*, ^ 7*^, 7*/ ^ 7*/ uniformly and £(7*,) ^ £(7xj)) as j 00. 
Hence, we see that 


i^(^.)-A7;)i < i^(r,)-A7:,)i + i^(7:,)-A7;)i 

< - >C(7x„)l ^0 as j ^ 00. 

Therefore, we have a subsequence of {rj}, still denoted by {vj}, that tends to a curve rjj 
uniformly as j —?> 00 with £(r(j) = £(7*^). Since 


II n - 7 ^# llco> - Ik# - rj llco + II Tj - 7*. Ilco - II 7*. - 7*^ |lco> 62, 

there are two minimizing curves, which is a contradiction. 

( 2 ) For the minimizing random walk 7 = 7(xn+i,0;^*) for k+i; have with ( 1 ) of 
Proposition 13.21 


•^n+1 




l'<k <0 


Ax 


+e( 7 ^, 4 -i)}At+ u( 7 '',b') + h{c){-ti>) 
^M(-;e)[ Y1 +^(k(7),4-i)}At 


l'<k <0 


+u(V'(7),b/) + h(c)(-£/) + e4(A) 


= E,,(. 




{L(‘=)( 7 a( 7 )(s), s, 7 a( 7 )'(s)) + ^( 7 a( 7 )(s), s)}ds 




+ h{c){-ti>) + e5(A) 


= E,,(.. 




E{va{i)) +e5(A). 


On the other hand, k_,_i = v{xn+i, 0) = £(7k J- Hence, we have 


( 4 . 16 ) 


0 < £, 




E{va{i)) - >C(7k+i) 


< 66(A). 


Here, e6(A) > 0 are independent of Xn+i- Set 

0 + := {7 G O I £( 7 a( 7 )) - ^i.lL+i) A e6(A)^}. 

Then, by (I 4 . 16 p . we have prohiVt'^) < e6(A)^/^. It follows from ( 1 ) of this lemma that 
there exists £o(^) > 0 for which, if |A| is small enough to realize e6(A)^/^ < eQ{u), we 
have II 7 a( 7 ) ~ 7 k+i llcoA ^ for all 7 G r 2 \r 2 ’''. Therefore, we obtain r 2 \r 2 + C hlAl^n+i) 
and 1 - e6(A)i/2 < pro&( 0 ^(x„+i)). 

Similarly, for the minimizing random walk 7 = 7(a:„+i, 0 ; for tcki’ have 
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with fl4.9p . fl4.12p and (3) of Proposition 13.1[ 


w. 


n+1 


V<k<0 
+ /ia(c)(— t/') 

l'<k<0 

+h(c){—tii) + e 7 (A) 


,0 


>C(r/A(7)) +e8(A), 


\w^+i - w„+i| = \w^+i - < ei(A). 

Therefore we have an estimate similar to fl4.16l) . and we may follow the same way as the 
above. □ 

Now, we compare fl4.8p and fl4.1ip . Take n so that is equal to 

max 

U^iS/3)\k=o 

Let T G N be such that qT > -^ log 4, let /' := —qT ■ 2K and let 7^^^^ : [—qT, 0] ^ M be 
the minimizing curve for n(x„+i,0). Then, fl4.15p implies 


hLJ-<lT)-j:{-qT)\ < -|7:„,,(0) -7r(0)|. 


Let be the minimizer for Then, it follows from fl4.8p and fl4.1ip that 

w°+i - w°+i < {/(7^ 4-i) - ^(7^ 4-i) + ^^(7^ 4-i) 


l'<k<0 


-£(7^4_l)}A^ + + (hA(c) - h{c))qT 


Ax 


l-qT 


7(7*') 

^(7L+i(s),s)c?s 


(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 


Ax 
Ax 


t-Tl+l 

>0 

'-J-qT 


MVA{l)(,s),s)ds 


Ax 


l-qT 


74(7x„+i(s),s)ds 


+ ^^M-;e)[ A(7^4-l)A^- / 7l(7A(7)(s),s)ds 

-^/7(';e)[ e(7^4-l)A^ 

l'<k<0 

+ -^/ 7 (-g*) (E {/(7^, 4 -i) - ^(7^, 4-i)}At 


+ -^/7(-g*) 


l'<k<0 

t V V 

' m(7‘ ) m(7‘ ) 

,0 


)-(m 


n+1 t^n+l) 


+ (w„+i - Vn+i) + (hA(c) - h{c))qT. 


The terms in (I4.18p - (l4.22p are denoted by i?i,...,i? 5 , respectively. Note that A{x,t) 
is uniformly continuous. By (2) of Lemma 14.31 and continuity of A, we obtain Ri < 
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eii(A)Ax. It follows from ( 1 ) of Proposition l3A] that 


prob{'y e I — 7^1 > Aa;4} < bsAx^, 


and hence we obtain i ?2 < e 9 (A)Aa:. Since e = o(Aa:), we have R 3 < eio(A)Ax. Since 
the minimizer (* is given as = Ilp(xm, tk, c + and satishes the Lipschitz 

condition in (2) of Proposition 13.21 Chebychev’s inequality yields for each a > 0, 


/\ ' 7 ' 

(4.23) probiy G hi I — 7^1 > «} < —^ for all /' < fc < 0. 

Since the averaged path 7 a of the minimizing random walk is C^-close to and 

tends to y*, the points (7 a(5),s) always stay in f/*(4(5/9) for small A. Hence, we 
have {(x, s) G I lx — 7 a(s)| < 5/2} c UH5). With the fact /(y^, tk-i)—g(y^, tfc-i) = 0 
on W{6), (Km and (023]) with a = 5/2, we obtain 


-R4 A ei2(A) 


Ax 

WW' 


Due to the choice of n, we have ~ ('^n+i ~ '^n+i) < 0 for 7 such that 7 *' 

is contained in the <5/3-neighborhood of y*{—qT). We see that the <5/6-neighborhood of 
7 a is contained in the (5/3-neighborhood of 7 *(—gT). In fact, it follows from fl4.17p 
that we have gP)— 7 /(—gP)| < <5/12, and hence we have | 7 a(— gT) — 7 *(—gT)| < 

1.5(5/12 for small A. By fl4.9p and fl4.23p with a = 5/6, we obtain 


R 5 A ei3(^) 


Ax 

j5j^' 


It follows from (I4.15p and the relation Vxx = Cxxi'^tt = Cn that 



-4(77,W,<i)*-rr‘ 


< v(S), 


where i>{5) comes from the modulus of the continuity of A. Therefore we obtain 

(4.24) (h(c) - Ai(c))«T < TP A + + e,,(A)Ax. 

We remark that if we do not introduce tCm+ii or if we do not choose the above n, it is 
not easy to see that P* is the principle term in (I4.24p . 

Next, we will obtain a converse estimate with the minimizer for where we 

take n so that — n°+i is equal to 

min 

t/7<5/3)U=o 
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It follows from fl4.8p and fl4.1ip that 


l'<k<0 

-e( 7 ^ 4 _l)}A^ + - /(y')) + (hAic) - h{c))qT 




2A J-qT 
Ax 

Ax 


'-J-qT 


Ax 


A{r]A{-f){s),s)ds - ^ 


qT 


l'<k<0 


l-qT 


A{T]A{-f){s),s)ds 


l'<k<0 

[E {/( 7 ^ tk-i) - ^/( 7 ^ 4 -i)}At 

V<k<Q 


+E^(,e) - / (y')) - (w^n+1 - ^^n+l) 

+ + (hA(c) - h{c))qT. 


Since the averaged path 7 ^ of the minimizing random walk is C'^-close to 
Jxn+i ts^ds to 7 *, the points ( 7 a(s),s) always stay in f/*(4(5/9) for small A. Hence, by 
(A4), the minimizer which is given by = Hp{xm,tk,c + DxW^_^_^), satishes the 

Lipschitz condition in ( 2 ) of Proposition 13.21 Therefore, with the estimate = 0(Ax), 
we can repeat the same argument as the above for the lower bound of each term. In 
particular, due to the choice of the n, we have “'^^( 7 *')) ~ ('^n+i ~'^n+i) > 0 for 7 

such that 7^' is contained in the (5/3-neighborhood of 7*(— gT), and the 5/6-neighborhood 
of 7 a (— gP) is contained in the 5/3-neighborhood of 7*(— gT). Thus we obtain 


(4.25) (h(c) - hA(c))qT > 7T*^ - “ ei 5 (A)Ax. 

In this way, we also obtain for 7 *, 

(4.26) (h(c) - h^{c))qT < + ei 5 (A)Aa:, 

(4.27) (h(c) - hA{c))qT > “ ei 5 (A)Aa:. 

If T* 7 ^ P and if we take small 5 > 0, there is a contradiction in fl4.24p . fl4.25p . fl4.26p and 

fl4.27p as A —)■ 0. This ends the proof of Proposition 14.21 □ 


Currently, we do not have any appropriate estimate of without (A4) for the rate 
of the law of large numbers in the minimizing random walk for A possible way to 

remove (A4) would be to sharpen Chebychev’s inequality. However, this seems not easy, 
because our random walk is not the sum of i.i.d. random variables. 


22 



























Next we show that, for certain values of c, the limit of any convergent subsequence 
has transition points at 7 **. For this purpose, we introduce the following symbols: 
Let x'l (resp. x'^) G {pr’ 7 *( 0 ),... be the nearest point to pr'j*,{k) on the 

right (resp. left) for fc = 0 ,..., q—1 (if 7 **( 0 ) = 0 , Xq is the point of {pr’ 7 *( 0 ),... ,pr'y*{q— 
that is nearest to 1). Let (reps. S^) be the segment between pr'^*„{k) and x^ 
of the projected unstable (resp. stable) manifold at f = 0 of 7 **. Set 

:= [area surrounded by and 5^]. 

0<k<q 

It follows from Subsection 3.3 that, for each c E J~ := [cq, ci — d~], there exists such 
that graph(c + v^'')\t=o 3 Uo<fc<gW^- Similarly, for each c E J~^ := [cq + (i+,Ci], there 
exists such that graph(c + hi'^^)|t=o ^ Uo<fe<gW^. Since there exists such that 
graph(c + hi''^)|t=o Uo<fc<g(W^ LI W^), we see that J~ fl J+ 7 ^ 0. 

Proposition 4.4. For each c E J~ (1 the limit of any convergent subsequence 
has transition points at 7 **. 

Proof. Let c G J” fl J"*". Suppose that there exists a convergent subsequence (ua = 
whose limit v does not have transition points at 7 **. The following is our strategy: Since 

V has transition points at some 7 * {i 7 ^ i*), we have the lower estimate fl4.25p . Since 

V does not have transition points at 7 **, the upper estimate fl4.24p with i = i* is not 
directly available. However, due to the choice of c and uniqueness of the value h(c), we 
have the Z^-periodic viscosity solution v = - 0 ^) of fll. 2 p such that '&( 7 **( 0 ), 0 ) = 0 and 

graph(c + h,j,)|4=o L U UM+). 

0<k<q 

Then, we can compare and v with modihcation of the argument in Proposition 14.21 
obtaining the upper estimate fl4.24p with i = i*. (A3) yields contradiction as A —)■ 0 . 

Adding constants if necessary, we have hA( 7 **( 0 ), 0 ) = h( 7 **( 0 ), 0 ) = 0. The entropy 
condition of Vx implies the following three cases: graph(c + Vx)\t=o contains 

(i) Sjf and a part of for A: = 0 ,..., g — 1 , 

(ii) a part of and for /c = 0 ,..., g — 1 , 

(iii) and S'f for A: = 0,..., g — 1. 

Take small 5 > 0 so that 

(i) ^ graph(c + h^) 3 prWL(<5)|{(x,t)gR2 | 77 (t)<a;< 77 (t)+ 5 }’ 

(ii) ^ graph(c + h,,) 3 

(iii) ^ 5 < min{mintgK |7r*-i(A) - 7** (f)I,minteK 17^*^ “ 7r*+iWI}- 

Following the proof of Proposition 14.21 we dehne a Z^-periodic C^-function v with h, 
instead of h, so that 

(i) ^ u = h on 1/(5) := {{x,t) E | 7**_i(A) <x< 7 **(A) + 5}, 

(ii) ^ V = V on V{6) := {{x,t) G | 7 **(A) - S <x < 7 **+i(A)}, 

(iii) ^ n = h on 1/(5) := {(x,A) G | 7 **_i(A) < x < 7 **_,_i(A)} (independent of 5). 
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Let denote the difference solution that yields va- We also define from 

and V with V{6) := V{6) fl Qodd as 

( h^+i on V (5) 

:= < v^+i + h^*(fc)+i - <*(A :)+1 for m*(k) < m, 

[ ^m+i + ^L(k )+1 - ^L(k )+1 for m < m,(k), 


where {xm»{k)+i,tk) and {xm^,{k)+i,tk) stand for the right end and left end of V{6)\t=t^, 
respectively. Then we may have g and /, as well as the difference equations for r’m+i 
v{xm+i,tk) and ^^+ 1 ; fo the same way as that of the proof of Proposition 14.21 

Since Vx = Vx on the boundary dV{S), we have for the ei 6 (A)-neighborhood of dV{6), 
denoted by B{A), 

sup \vx -Vx\< €n{A). 

B{A) 

Since (hA)x W uniformly on V{6) including its boundary as A —)■ 0 , it follows from the 
continuity of Vx on dV (<5) and the relation |(hA)a; — r’a;| < |(hA)x — + \vx — Vx\ + \vx — Vx\ 

that we have 

sup |(hA)a; “ Vx\ < ei 8 (A). 

B(A) 

Hence, fl4.13p and fl4.14p imply that fl4.12p holds in all the cases (i)-(iii). Note that we 
fail to have fl4.12p . if we simply follow the proof of Proposition 14.21 with t/® (5) in the 
definition of We now do not have the estimate < ei 9 (A) on Qodd\k=Q- 

However, due to the choice of v and the normalization of h, h, va at 7 **( 0 ), it follows from 
the relation v{x, 0 ) — v{x, 0 ) = ^^•^{vx{y, 0 ) — Vx{y, ^))dy that we have 


(i) v{x^ 0 ) — ' 0 (x, 0 ) 

(ii) ^ h(x, 0 ) — v{x, 0 ) 


0 for 7**(0) < X < 7**(0) + 6, 
negative for 7 **_i( 0 ) < x < 7 **( 0 ), 

0 for 7**(0) — 5 < X < 7**(0), 
negative for 7 **( 0 ) < x < 7**^1 ( 0 ), 


(iii) v{x, 0 ) — v{x, 0 ) 


0 forx = 7 **( 0 ), 

negative for 7 **_i( 0 ) < x < 7 **_,_i( 0 ) and x 7 ^ 7 **( 0 ). 


Since va v uniformly as A —)■ 0 and 


= V^^i-v(Xm+l,0) 

< ei 9 (A) +v(Xm+l,0) -v(Xm+l,0) onV((5), 

i,P - iP = tP - iP 

a^m+l ^m+1 ^m*(0)+l ^m*(0)+l 

< ei 9 (A) + t;(a;m*(o)+i, 0 ) - v{xm*{o)+i, 0 ) for ni*{k) < m, 

i,P - iP - tP - iP 

^m+l ^m+1 ~ ®^m*(0)+l ®^m*(0)+l 

< ei 9 (A)+n(xm 40 )+i, 0 )-f)(xm 40 )+i, 0 ) form<m*(fc), 

we have in all the cases (i)-(iii), 

W^m+l - < ei9(A) on Qodd\k=0- 

Take n so that is equal to 

max 

(7»*(5/3)U=o 
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Here, is not Note that 

-620 (A) < 


because 

^n+i-'^n+i = hA(x„+i,0)-h(a;„+i,0)+h(x„+i,0)-t)(a;„+i,0) 

> V^{Xm,{'y*,{0))^ 0) — 'u(Xm,(7** (0))) 0) + v(Xm(^*^ (0)) ^ 0) — v(Xm('j*^ (0)) y 0), 

wa — h —)■ 0 on h"**((5/3)|fe=o and h( 7 **( 0 ), 0 ) — '&( 7 **( 0 ), 0 ) = 0. Then, we can estimate 
tc°+i — from the above with the minimizer for in the same manner as that 
of the proof of Proposition 14.21 obtaining 

(4.28) (h(c) - hA{c))qT < TT"^*^ + e 22 (A)Ax. 

On the other hand, by the argument in the proof of Proposition 14.21 we have for 7 * 
{i ^ i*) at which v has transition points, 

(ft(c) - fti(c)),T > TV^ - v{S)^ - ‘-^Ax - £24(A)Ax. 

Thus, by (A3), we reach a contradiction as A ^ 0. □ 


We extend Proposition 14.41 

Proposition 4.5. For each c G (co,ci), the limit of any convergent subsequence 
has transition points at 7 **. 

Proof. Let c* := max J“ fl J'^. Then, Proposition 14.41 implies that the limit of any 
convergent subsequence has transition points at 7 **. By (A3) and Proposition 

14.21 has transition points only at 7 **. It follows from Subsection 3.3 and the choice 
of c*that we may have the Z^-periodic viscosity solution such that graph(c* + 
hi'" ^)| 4 =o F> Uo<k<q{h{k graph(c* + hi'" '^)\t=o \ ^o<k<qh{k is a part of the upper 

separatrix at t = 0. It is clear that has transition points only at 7 **. By Proposition 
an and must coincide (up to constants) with the Peierls barrier. Therefore 
we see that 

graph(c* + hi'" i)|t=o =graph(c* + hi'" '')\t=o- 
Suppose that there exists c G (c*,ci) for which (*) does not hold: 

(*) The limit of any convergent subsequence has transition points at 7 **. 

Then, we have a convergent subsequence whose limit does not have transition 

points at 7 *,. It follows from (2) of Proposition 13.11 that we have c* + (h^ < c + 

(hA^)x a.e., obtaining c* + Vx ^|t=o < c + Vx^\t=o a.e. Hence, if graph(c + does 

not contain , then must have transition points at 7 **. Therefore, we see that 
graph(c + hi^^)|f=o D Uo<k<qSk . However, this means that graph(c + hi‘^^)|f=o coincides 
with the whole upper separatrix at t = 0, which is only allowed for c = ci. Therefore, 
(*) holds for all c G (c*, Ci). 

We have the lower extension with c* := min J~ fl in the same manner. □ 
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Proof of Theorem \2. li By Proposition 14.51 the limit of any convergent subsequence 
must have transition points at 7**. By (A3) and Proposition 14.21 has transition 
points only at 7**. By Proposition 14.11 such is unique up to constants and co¬ 
incides with hp‘^^(7*»(0), 0; •, •). Therefore, with the normalization, we conclude that 
h.p^^(7**(0), 0; •, •) as A —)■ 0 in the whole sequence. By (5) of Proposition 13.11 
this implies that Oj■))a: as A ^ 0 in the whole sequence. □ 
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